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Abstract: 

The associate production of Higgs bosons with W or Z bosons, known as Higgs-strahlung, 
is an important search channel for Higgs bosons at the hadron colliders Tevatron and LHC for 
low Higgs-boson masses. We refine a previous calculation of next-to-leading-order electroweak 
corrections (and recalculate the QCD corrections) upon including the leptonic decay of the W/Z 
bosons, thereby keeping the fully differential information of the 2-lepton -|- Higgs final state. 
The gauge invariance of the W/Z-resonance treatment is ensured by the use of the complex- 
mass scheme. The electroweak corrections, which are at the level of —(5—10)% for total cross 
sections, further increase in size with increasing transverse momenta in differential cross 
sections. For instance, for pt,h ^ 200 GeV, which is the interesting range at the LHC, the 
electroweak corrections to WH production reach about —14% for Mh = 120 GeV. The described 
corrections are implemented in the Hawk Monte Carlo program, which was initially designed 
for the vector-boson-fusion channel, and are discussed for various distributions in the production 
channels pp/pp — )■ H -|- lui / 1~ l^ / vii^i + X. 
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1 Introduction 



The Higgs-strahlung processes at hadron colliders, pp/pp — )• H + W/Z — )• H + Ivi / / uiui + X, 
represent important search channels for a Standard Model (SM) Higgs boson in the low-mass 
range, which is not only favoured by electroweak precision tests, but also by recent results from 
Higgs searches at the LHC [Ij. At the Tevatron the sensitivity to a light SM Higgs boson almost 
exclusively results from these channels, and the analysis of the full Tevatron data set might be able 
to exclude values of the Higgs-boson mass Mh close to the LEP exclusion limit of Afn = 114 GeV 
beyond the 95% confidence level. At the LHC this low-mass range is explored mainly via inclusive 
Higgs production employing the H — )• 77 decay channel, but Higgs-strahlung slightly improves 
the sensitivity and provides an interesting additional independent search channel with the Higgs 
boson decaying via H — )• bb. The challenge in all these analyses is background control, both in 
cross-section normalization and differential distributions. This is particularly pronounced in the 
boosted H + W/Z analysis at the LHC where the H — )• bb decay should be identified from the 
b-quark substructure within a highly boosted "fat jet" [2|. 

The mentioned analyses call for precise theoretical predictions with proper estimates of the 
corresponding theoretical and parametric uncertainties. The next-to-leading-order (NLO) QCD 
corrections, whose structure is identical to the respective corrections to the Drell-Yan process, 
were calculated ^ a long time ago and are available, e.g., in the public programs VV2H ^ and 
MCFM [5]. The most important part of the next-to-next-to-leading-order (NNLO) QCD correc- 
tions is similar to the Drell-Yan process as well and was calculated in Ref. [6], using the results 
of Ref. [7] on the total Drell-Yan cross section. The production of ZH final states additionally 
receives contributions from gluon fusion via top-quark loops, which was also included in the 
calculation of Ref. [6] . The remaining NNLO QCD corrections that are not of Drell-Yan type, 
which involve Higgs couplings to top-quark loops, have been shown to be at the 1—2% level [8]. 
The calculation of NLO corrections, which are of the order of +30% in QCD, was completed by 
the NLO electroweak (EW) corrections ^J, which turned out to be about —(5—10)%, and a com- 
bined prediction for total cross sections at the Tevatron and the LHC was presented in Ref. |10j . 
An update of the cross-section prediction, together with a thorough estimate of uncertainties, 
was presented in the first report [11] of the LHC Higgs Cross Section Working Group. For the 
LHC with a centre-of-mass (CM) energy of 7(14) TeV the QCD scale uncertainties were assessed 
to be about 1% and 1—2(3—4)% for WH and ZH production, respectively, while uncertainties of 
the parton distribution functions (PDFs) turn out to be about 3—4%. 

These results, however, hold only for total cross sections and have to be generalized to 
distributions for the experimental analyses, in particular if one employs boosted Higgs bosons. 
An important step towards differential distributions has been made recently at NNLO QCD in 
Ref. [12], where the Drell-Yan-like corrections to WH production are calculated while keeping 
the full kinematical information of the final state. 

In this paper we present fully differential NLO EW predictions for the final states H + 
Ivi/l^l^ /viUi + X, i.e. including the leptonic decays of the W and Z bosons. We also recalculate 
the corresponding NLO QCD corrections. In contrast to the (N)NLO QCD calculations the 
EW loop corrections involve also irreducible corrections connecting production and decay of the 
W/Z bosons. In order to guarantee the gauge invariance of the NLO calculation in spite of the 
W/Z resonances we employ the complex-mass scheme, which was introduced for leading-order 
(LO) and NLO calculations in Ref. [13] and Ref. |14] . respectively. The use of complex W/Z 
masses at the same time cures the appearance of perturbative artifacts for Higgs-boson masses 
near the pair production thresholds at Mh = 2Mw, 2Mz, where the EW corrections to WH/ZH 
production [9] are singular. The calculation described in this paper will be part of the Monte 
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Figure 1: Feynman diagrams for the LO processes (|2.ip - ()2.4p . 



Carlo program Hawk, which was originally designed for the description of Higgs production via 
vector-boson fusion including NLO QCD and EW corrections [15j and is publically available [16] . 

The paper is organized as follows: In Section [2] we describe the structure of the underlying 
NLO calculation and the techniques used. Section [3] contains a detailed discussion of numerical 
results for the processes pp/pp — )• H + Zf;//~/+/z/;P; + X at the Tevatron and the LHC, the latter 
at CM energies of 7TeV and 14TeV. Finally, we conclude in Section HI 



2 Structure of the NLO calculation 
2.1 General setup 

At LO, the associate production of a Higgs boson H with a weak gauge boson V = W, Z can only 
proceed via quark-antiquark annihilation at hadron colliders. Treating the incoming quarks and 
the outgoing leptons as massless, the Higgs boson does not couple to the massless fermions, and 
there is only one LO diagram per channel, see Fig. [H In detail, the following partonic processes 
are considered, 

Ui dj HW+ ^ , (2.1) 

dj Ui HW- m-vi , (2.2) 
qi Qi ^ HZ ^ Hr/+ , (2.3) 
Qi qi ^ HZ ^ Ylvivi , (2.4) 

where qi denotes any light quark and Uj, di the up- and down- type quarks of the ith generation. 
The intermediate W/Z-boson resonances are described by complex W/Z-boson masses /iy via 
the replacement 

^ /x^ = - iMyFy, y = W, Z (2.5) 

in the V propagator as dictated by the complex-mass scheme (see below). Hence, all our results 
correspond to a fixed-width description of the Breit-Wigner resonance. Moreover, all related 
quantities, in particular the weak mixing angle, are formulated in terms of the complex mass 
parameters. 

The final-state leptons are treated as massless unless their small masses are used to regularize 
a collinear divergence. Concerning bremsstrahlung, we support the possibility that collinear 
photons may be completely separated from an outgoing charged lepton, because this situation 
is relevant for muons. More details on the treatment of such non-collinear-safe observables are 
described below. 

The light quarks are considered massless as well, in line with the parton-model requirements. 
This means that the quark mixing matrix only appears as global weight factor \Vij\'^ in the 
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processes (j2.ip and (j2.2p . Mixing with the third generation is neglected throughout. Owing to 
the smallness of the bb contribution to ZH production, which amounts to approximately 1(3)% 
at the 7(14) TeV LHC for an inclusive analysis, 0.5(1)% for a boosted analysis, and is negligible 
at the Tevatron, the bb annihilation channels are treated in LO only. 

To define the electromagnetic coupling constant a, we use the scheme, i.e. we derive a 
from the Fermi constant according to 
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In this scheme, the weak corrections to muon decay, Ar, are included in the charge renormal- 
ization constant (see e.g. Ref. [9j). As a consequence, the EW corrections are independent of 
logarithms of the light-quark masses. Moreover, this definition effectively resums the contribu- 
tions associated with the running of a from zero to the weak scale and absorbs some leading 
universal corrections oc G^m^ from the p parameter into the LO amplitude. For corrections due 
to collinear final-state radiation it would be more appropriate to use a(0) defined in the Thomson 
limit to describe the corresponding coupling. On the other hand, using everywhere is best 
suited to describe the larger genuine weak corrections which are present in the Higgs-strahlung 
processes. Thus, the optimal choice cannot be achieved in one particular input scheme, and the 
calculation could be refined beyond NLO. Among other things, higher-order effects from multi- 
photon emission should also be included at this level of precision which is beyond the scope of 
this work. However, this kind of sophistication will not be needed in practice, since the expected 
statistics does not allow for this level of accuracy in the experimental analysis either. 

The major part of the calculation could be carried over from the calculation of NLO QCD 
and EW corrections to Higgs + 2 jet production as described in Ref. [15]. In detail, we had to 
discard the vector-boson fusion contribution (t/u channels) and to keep only the s-channel part 
from Ref. [15] adjusting the quantum numbers of the final-state fermions, which are leptons for 
the Higgs-strahlung processes under consideration. 

All results are checked against a completely new and independent second calculation based on 
FeynArts 3.2 [17] and FormCalc 3.1 [18]. The translation of the amplitudes into the Weyl- 
van-der-Waerden formalism as presented in Ref. [19] is performed with the program Pole |20j . 
Pole also provides an interface to the multi-channel phase-space integrator Lusifer [21] which 
has been extended to use Vegas in order to optimize each phase-space mapping. The results of 
the two calculations are in mutual agreement. 



2.2 Virtual corrections 

The actual loop calculation, which was also the basis for the loop calculation [15] entering Hawk, 
goes back to the calculation [22] of NLO QCD and EW corrections to the four-body Higgs decays 
H —7- WW/ZZ —7- 4 fermions, where the same amplitudes appear in a crossed variant. We make 
use of those results, which were obtained in the traditional Feynman-diagrammatic approach and 
checked against the independent second calculation. 

The virtual corrections modify the partonic processes that are already present at LO; there 
are about 200 one-loop diagrams per tree diagram depending on the channel. At NLO these 
corrections are induced by self-energy, vertex, box (4-point), and pentagon (5-point) diagrams, 
shown schematically in Fig. [2l For the implementation of the finite W/Z widths the complex- 
mass scheme is employed, which was introduced in Ref. [13] for LO calculations and generalized 
to the one- loop level in Ref. [H]. In this approach the W- and Z-boson masses are consistently 
considered as complex quantities, defined as the locations of the propagator poles in the complex 
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Self-energy diagrams: 




Figure 2: Different classes of loop diagrams with up to 5 external legs corresponding to 
the virtual corrections to the LO processes (|2.ip - (|2.4p . Note that there are loop-induced 
vertex diagrams for the coupling of the Higgs boson to the external massless fermions 
which are absent at tree level. 



plane. This leads to complex couplings and, in particular, a complex weak mixing angle. The 
scheme fully respects all relations that follow from gauge invariance. 

The underlying calculation |22] of the H — )• 4/ amplitudes was performed both in the con- 
ventional 't Hooft-Feynman gauge and in the background-field formalism using the conventions 
of Refs. [23] and [24], respectively. The amplitudes were generated with FeynArts 1.0 [25] 
and further manipulated with in-house programs written in Mathematica, which automatically 
create building blocks in Fortran. The amplitudes are expressed in terms of standard matrix 
elements, which contain the fermionic spinor chains, and coefficients, which contain the tensor 
integrals, couplings, etc. The tensor integrals are recursively reduced to scalar master integrals 
at the numerical level. Tensor and scalar 5-point functions are directly expressed in terms of 
4-point integrals 126^127). Tensor 4-point and 3-point integrals are reduced to scalar integrals 
with the Passarino-Veltman algorithm [28] as long as no small Gram determinant appears in 
the reduction. If small Gram determinants occur, we expand the tensor coefficients about the 
limit of vanishing Gram determinants and possibly other kinematical determinants, as described 
in Ref. |27j in detail. For the evaluation of the scalar integrals with complex masses we use two 
independent in-house libraries based on the general results listed in Ref. |29j . 
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Figure 3: Feynman diagrams for real photonic bremsstrahlung corrections to the LO 
process (j2.ip . A photon cannot only be emitted from the external charged fermions, but 
also from the intermediate W bosons. In Feynman gauge, there are also diagrams involving 
the charged Goldstone boson 



2.3 Real corrections 

The EW real corrections to the partonic processes (I2.ip - (l2.4p consist of photon bremsstrahlung, 



with Feynman diagrams shown exemplarily for the process (j2.7p in Fig. [3l and of photon-induced 
processes, which involve the same amplitudes as the bremsstrahlung processes, but with the 
photon crossed into the initial state in all possible ways. The technical implementation of the 
real corrections again basically follows along the same lines as in the NLO calculation within 
Hawk for Higgs + 2 jet production ^15iil6j, employing dipole subtraction [30H32] . 

Due to the emission of soft photons the real corrections include soft singularities which are 
cancelled by the virtual corrections independently of the details of the event selection or recombi- 
nation procedure. The treatment of left-over collinear singularities related to splittings involving 
photons and fermions requires particular care: 

• Photonic initial-state radiation 

The collinear singularities from photonic initial-state radiation off quarks are absorbed by 
a redefinition of the PDFs similar to the standard factorization procedure in QCD, but 
employing the DIS scheme. At present, the PDF set MRSTQED2004 [33] is the only one 
that takes into account 0{a) corrections. However, applications with this PDF set and 
older estimate of 0{a) effects on PDFs show that the impact of these PDF corrections 
is at most at the per-cent level [34j. For the time being, it is, thus, advisable to use an 
up-to-date PDF set that ignores 0{a) corrections, but includes recent improvements on 
the QCD side, because the latter have an impact on PDFs much larger than 1%. 

• Photon-induced processes, 'y^qq*/qq* splitting 

In the photon-induced processes, the collinear singularities resulting from the splittings 



Ui dj HW+ ^ iiuil+-f , 

dj Ui — ^ HW^ — ^ }il~i?ij, 

qi ^ HZ ^ Hr /+7 , 

qi Qi ^ HZ Hi^ii^n , 



(2.7) 
(2.8) 
(2.9) 
(2.10) 
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7 ~^ QQ*/qQ* s-re removed in the process of PDF redefinitions as well. Moreover, the 
evaluation of the photon-induced channels requires a photon distribution function and its 
inclusion in the PDF evolution. Using the MRSTQED2004 PDF set for this task is certainly 
legitimate, since the contribution from initial-state photons is comparably small. 

• Collinear final-state radiation 

Next, we consider collinear final-state radiation off leptons. If collinear photons and charged 
leptons are recombined into a pseudo-particle (mimicking the start of hadronic or electro- 
magnetic showers) to form IR-safe observables, all the remaining singularities arising from 
collinear photon emission in the final state also cancel against the corresponding singular- 
ities in the virtual corrections. For muons in the final state, however, it is experimentally 
possible to separate collinear photons from the lepton, i.e. to observe so-called "bare" 
muons. Hence, the resulting cross sections are not collinear safe (i.e. the KLN theorem [35] 
does not apply), and the corresponding collinear singularities show up as logarithms of the 
small lepton (muon) mass. The lepton mass cuts off the collinear divergence in a physi- 
cally meaningful way. We employ the extension [32j of the subtraction formalism [301131]. 
which allows for calculating cross sections for bare leptons, i.e. cross sections defined with- 
out any photon recombination. Like in the standard subtraction formalism, it is sufficient 
to calculate the real-emission matrix elements for the partonic processes in the massless- 
fermion approximation. The main difference between the subtraction variants of Refs. [31] 
and [32] concerns the implementation of phase-space cuts. In the standard subtraction 
formalism [31] it is always assumed that the complete momentum of the lepton and a 
collinear photon is subject to cuts (as it would be the case after recombination), while 
the generalization [32] allows for non-collinear-safe cuts that resolve the distribution of the 
momenta in collinear photon-lepton configurations. This more general cut procedure in 
the non-collinear-safe case has to be carefully implemented both in the real-emission part 
and the corresponding subtraction terms, in order to ensure the numerical cancellation of 
singularities. Of course, the treatment of non-collinear-safe cuts leads to a modification 
of the re-added subtraction part as well. In the formulation of Ref. [32j this modification 
assumes the form of an additional (-l-)-distribution which contains the surviving mass sin- 
gularity. This (-l-)-distribution integrates to zero for collinear-safe observables so that the 
formalism reduces to the well-known standard subtraction formalism. For non-collinear- 
safe observables the additional logarithms of the lepton mass in the final result are, thus, 
isolated analytically. For a complete and detailed description of this more general subtrac- 
tion formalism we refer the reader to Ref. [32] . 

• Photon-induced processes, 7 — ?■ II* /II* splitting 

Finally, there are also collinear configurations connected with the splittings 7 — t- II* /II* if 
charged leptons are allowed to escape into the beam pipe. For instance, in our numerical 
analysis in Section [3l we also provide predictions for a contribution to the uivi chan- 
nels which stems from WH production, where the charged lepton is not identified in the 
detector. While the EW corrections are conceptually not different for the quark-induced 
processes, the photon-induced processes have to be treated carefully, in particular the con- 
tribution where the incoming photon is splitting into a charged lepton pair, which leads to 
an enhancement by a logarithm in the lepton mass. The subtraction formalism to extract 
these logarithms analytically is available [32j and proceeds as for the photonic splitting 
into quarks. However, in contrast to the splitting into quarks the mass logarithm is not 
absorbed into the definition of the PDF (there is no lepton PDF), but is part of the re- 
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suit. For the numerics in Section [3l we use the mass of the electron in this logarithmically 
enhanced term. 

3 Numerical results 

3.1 Input parameters and setup 

For the numerical evaluation we adopt the input parameters used by the LHC Higgs Cross Section 
Working Group in Ref. 

0.1202, 
2.0887 GeV, 

2.4952 GeV, MH = 120GeV, 

105.658367MeV, mt = 172.5 GeV, 

\Vcd\ = Vl-I^csl^ (3.1) 

which essentially follow Ref. [36]. Note that we use the measured width F^^ for the Z boson, 
but the calculated on-shell width F^^ for the W boson as input. The CKM matrix is included 
via global factors in the partonic cross sections for the different possible quark flavours. The 
small mixing between the first two and the third generation is neglected. Within loops the CKM 
matrix is set to unity, because its effect is negligible there. 

Using the complex-mass scheme |14) . we employ a fixed width in the resonant W- and Z-boson 
propagators in contrast to the approach used at LEP and Tevatron to fit the W and Z resonances, 
where running widths are taken. Therefore, we have to convert the "on-shell" (OS) values of 
My^ and Ty^ {V = W, Z), resulting from LEP and Tevatron, to the "pole values" denoted by 
Mv and Ty- The relation between the two sets of values is given by [37] 

Mv = MOS/v/l + (rOS/MOS)2, rv = T^^ / ./iUT^/MO^ , (3.2) 

leading to 

Mw = 80.370 . . . GeV, Fw = 2.0880 . . . GeV, 

Mz = 91.153 . . . GeV, Fz = 2.4943 . . . GeV. (3.3) 

We make use of these mass and width parameters in the numerics discussed below, although the 
difference between using My or My^ would be hardly visible. 

As explained in Section [2.H we adopt the scheme, where the electromagnetic coupling a 
is set to ac^ (12. 6p . In this scheme the electric-charge renormalization constant does not contain 
logarithms of the light-fermion masses, in contrast to the a(0) scheme, so that the results become 
practically independent of the light-quark masses. 

The MSTW2008NLO/LO PDF set [38] is used throughout implying the values of as(Mz) 
stated in (13. ip at NLO. We use standard running of the strong coupling constant as provided 
by the LHAPDF implementation of the PDF sets [39]. Only the photon-induced processes are 
evaluated with the MRSTQED2004 set of PDFs [33] • The QCD and QED factorization scales 
as well as the renormalization scale are always identified and set to 

H = Mv + Mu. (3.4) 





= 1.16637 X lO^^GeV-^, 


as(Mz) NLO 




= 80.398 GeV, 


rOS 




= 91.1876 GeV, 


rOS 
^ z 


me 


= 0.510998910 MeV, 




\Vud\ 


= iKsl = 0.974, 


\Vns\ 
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For a study of the scale and PDF uncertainties in Higgs-strahlung processes, we refer the reader 
to the report of the LHC Higgs Cross Section Working Group [TT] for total cross sections and to 
the upcoming second report ^40j on differential distributions. 

In the following, we present results for the Tevatron at a CM energy of y/s = 1.96 TeV and 
for the LHC with = 7 TeV and 14 TeV. 

3.2 Phase-space cuts and event selection 

Since WH and ZH production with leptonically decaying W and Z bosons do not involve outgoing 
QCD partons in LO, there is at most one QCD parton in the final state in our fixed-order NLO 
predictions, so that no jet algorithm has to be applied. 

As explained above, we alternatively apply two versions of handling photons that can become 
collinear to outgoing charged leptons. The first option is to assume a perfect isolation between 
charged leptons and photons, an assumption that is at least approximately fulfilled for (bare) 
muons. The second option performs a recombination of photons and nearly collinear charged 
leptons and, thus, mimics the inclusive treatment of electrons within electromagnetic showers 
in the detector. Specifically, a photon 7 and a lepton / are recombined for R^i < 0.1, where 

Rryi = ^ (yi — y^)^ + (pf^ is the usual separation variable in the y— i;^-plane with y denoting the 
rapidity and (pi^ the angle between / and 7 in the plane perpendicular to the beams. If / and 
7 are recombined, we simply add their four-momenta and treat the resulting object as a quasi- 
lepton. If more than one charged lepton is present in the final state, the possible recombination 
is performed with the lepton delivering the smaller value of R^i. 

After applying the recombination procedure we impose the following cuts on the charged 
leptons, 

I^TJ>20GeV, \yi\< 2.5, (3.5) 

where pT,i is the transverse momentum of lepton I. For channels with at least one neutrino in 
the final state we require a missing transverse momentum 

> 25 GeV , (3.6) 

which is equal to the total transverse momentum of the neutrinos in the events. While the 
two-body decay of the Higgs boson is implemented in Hawk, in this work we do not employ 
identification cuts on the decay products of the Higgs boson. 

For the LHC, we also discuss the impact of the optional additional cuts 

PT,H > 200 GeV, Pt,w/z > 190 GeV (3.7) 

on the transverse momentum of the Higgs and the weak gauge boson, respectively. The corre- 
sponding selection of events with boosted Higgs bosons is improving the signal-to-background 
ratio in the context of employing the measurement of the jet substructure in H — )• bb decays 
leading to a single fat jet. The need for background suppression calls for (almost) identical cuts 
on the transverse momentum of the vector bosons and the Higgs boson [4T]. However, symmetric 
cuts induce large radiative corrections in fixed-order calculations in the corresponding distri- 
butions near the cut. Since the Higgs boson and the vector boson are back-to-back at LO, any 
initial-state radiation will either decrease pt,h or Pt,w/Z! the event may not pass the cut 
anymore. Hence, the differential cross section near the cut is sensitive to almost collinear and/or 
rather soft initial-state radiation. By choosing the above (slightly asymmetric) cuts this large 
sensitivity to higher-order corrections can be removed for the important pT,H-distribution. Of 
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course, since the LO distribution for Pt,w/z is vanishing for Pt,w/z < 200 GeV due to the pt,h 
cut, the higher-order corrections to the Pt,w/z distributions are stih large in this region. 

The channel with a Higgs boson and only missing also receives a contribution from WH 
production where the charged lepton is not identified, which we label Hf^/T/. In our setup, a 
lepton is not identified if it fails to pass either the pt,« or the yi cut in (|3.5p . Following our 
generic setup, we still consider only the neutrino momentum to calculate ^t. Since in this case 

and pt,w ars indistinguishable, we also discard the momentum of the charged lepton when 
calculating the transverse momentum of the vector boson, i.e. we set pt,v = in this particular 
case. 

3.3 Results for the LHC and the Tevatron 

In this section, we present numerical results for (differential) cross sections and various corrections 
to the Higgs-production channels with one charged lepton and missing transverse momentum 
(HZ+f/ and Hl~!7;), with two charged leptons (H/+Z^), and with missing momentum only in the 
final state. For the latter, we distinguish the Z-mediated process {HviVi) and the W- mediated 
(Kh'i/T'i) process, where the charged lepton is not measured in the detector. Of course, the 
two sources are not distinguishable experimentally. However, since the contribution of the two 
sources varies considerably with the applied identification cuts, we state the results separately. 
For the Hi'i/Vi channel we simply invert the identification cuts (j3.5p for the charged lepton 
(see discussion at the end of Section 13. 2p while all other cuts and the recombination procedure is 
applied without changes (see Section [2^ for subtleties concerning the photon- induced processes). 
Thus, the resulting cross section may serve as an upper bound on the contribution to the channel 
with a Higgs boson and only missing transverse momentum in the experimental analysis. 

Concerning the EW corrections, we distinguish bare muons and final states with lepton- 
photon recombination, as introduced in Section 13.21 The corresponding corrections are labelled 
^eW ^^'^ ^EW' respectively. While the latter is independent of the specific charged lepton in 
the final state, contains fermion-mass logarithms and depends on the lepton mass which is 
taken to be the muon mass. The contributions due to photon-induced processes are also given 
in terms of their relative size with respect to the LO cross section and are labelled 5^. 

All cross sections are given for a specific leptonic state (e.g. e+e", /x"'"/U~, or v^u^). While 
the bare corrections are only applicable for the muon final state, the results with lepton-photon 
recombination are equally valid for electrons and muons. For the process HuiVi, there is of 
course no recombination to be applied, i.e. 5|f^ = S^^/^ and the results are valid for all three 
generations. Hence, the total ZH cross section with invisible Z decay can be trivially obtained 
by multiplying the given result by three. On the other hand, the photon-induced contribution to 
Hi'l/T'l depends on the fermion generation and we show numerical results for the electron channel. 
The corrections in the muon and tau channels are smaller since the corresponding logarithms of 
the fermion masses are smaller. 

In a single run Hawk employs a single set of PDFs and produces the following NLO predic- 
tion, 

f^HAWK = <T0 X (1 + JqcD + ^7 + '^EW) • (3.8) 

Note that the LO part (Tq is calculated with NLO PDFs in that case, i.e. it is not equal to the 
proper LO cross section cr^^, which is to be calculated with LO PDFs and requires a separate run 
of Hawk. Likewise, (1 + (5qcd) is not equal to the standard QCD K-factor Kqqd = '''qcd/^^*^- 
In the spirit of factorizing QCD and EW corrections, which is certainly true for corrections that 
are dominated by soft or collinear gluons, the fixed-NLO prediction (|3.8p can be improved by 

^fact° = ^QCD X (1 + Sew) + 6„ a^hu = (1 + '^qcd) , (3.9) 
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and we show predictions for the electron case, i.e. Sew = ^ew ™ numerical results below. The 
NLO EW correction factor 5ew is rather insensitive to changes in the PDFs, but involves a small 
dependence on the QED factorization scale. Recall that the presently missing 0{a) corrections 
to the PDFs induce an uncertainty at the level ^ 1%, i.e. for full 0{a) precision QED-improved 
PDFs are needed. To fix full 0{aas) precision beyond this factorization ansatz, a full two-loop 
calculation would be required, which would certainly be overkill for the data analyses at the 
Tevatron and LHC. Equation (j3.9p also offers a simple possibility to further improve the NLO 
result to include the NNLO QCD predictions by other authors [6l[l2] upon replacing ctqq§ by 
'^QCD*^- production this is done in Ref. [30], where our results for ^ew are used. 

In Figure m we show our best predictions cr^^ for the cross section and the EW corrections 
for the five Higgs-production channels under consideration. We show the results for the inclusive 
analysis at the Tevatron as well as for the boosted analysis at the LHC running at a CM energy of 
7TeV. As expected, the EW corrections become more and more negative with rising CM energy 
of typical events and, thus, are larger in absolute terms for the boosted analysis. Their size 
depends on the channel and exceeds —15% for the boosted analysis. As a function of the Higgs- 
boson mass, the corrections show small dips around Mh = 2Mw and Mh = 2Mz corresponding 
to threshold singularities of vector-boson pairs inside loops which are automatically regularized 
by the complex-mass scheme in our calculation. Earlier predictions in Ref. [9] have been singular 
at these thresholds. Away from these thresholds, the corrections for the inclusive analysis differ, 
as expected, only at the level of 1—2% from the predictions for total cross sections for on-shell 
gauge bosons shown in Ref. [9], mainly due to the effect of final-state radiation in the presence 
of lepton identification cuts. For the inclusive analysis, the cross section of the Hz^^/f^ channel 
contributes a sizable fraction to the Higgs plus missing transverse momentum cross section. For 
the boosted analysis, where the vector boson has to possess a large transverse momentum, it is 
more unlikely for the charged lepton not to be detected and the contribution due to a missed 
charged lepton is negligible (below 0.1 fb at the 7TeV LHC). 

In Figure m we show the EW corrections for Hl'^i'i production at the Tevatron and the LHC 
at 7 and 14TeV. The differences between the corrections at the colliders are rather small, at 
the 1% level. As already indicated above, the EW corrections are insensitive to the choice of 
the PDF set. Note that there is only a small difference of about 1% or 2% between the cases of 
bare or recombined leptons, where the larger value of 2% refers to the boosted-Higgs analyses. 
This small difference reflects the fact that the major part of the EW correction is not due to 
final-state radiation, but due to non-universal weak corrections. 

In Figure El we show the contribution from photon-induced processes to the W-mediated 
channels at the LHC. The contribution to the channels, where the charged lepton is identified, is 
at the 1—2% level for the boosted analysis and mainly results from diagrams where the initial- 
state photon couples to a f-channel W boson. Accordingly, the corrections are almost zero 
for the Z-mediated channels (not shown). There are larger photon- induced corrections for WH 
production, rising to 8%, to the inclusive cross section at the LHC for large Higgs masses around 
200 GeV. Moreover, the photon-induced corrections to the Hi'i/T'i channel, which are additionally 
enhanced by logarithms of the charged lepton mass, for boosted Higgs bosons are larger than 
10% (for the photon splitting into an e''~e~ pair) of the corresponding LO Hz^^/F^ cross section 
and even more sizable for the inclusive setup, in particular for high Higgs masses. Whether these 
large corrections become phenomenologically relevant at the LHC depends on the details of the 
analysis. At the Tevatron the corresponding corrections do not exceed 3%. 

To facilitate a possible comparison of our results with future results of other groups, in 
Tables [T] and [21 we show the cross sections and corrections for Mh = 120 GeV. 
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Figure 4: The cross section crf^^^ (left) and NLO EW corrections (right) for the different Higgs- 
strahlung processes with basic cuts at the Tevatron as well as for boosted Higgs bosons at the 
7TeV LHC function of the Higgs-boson mass. 
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channel 


m+i^i + X 


m-i)i + X 


m+i- + X 


Ruiiyi + X 


+ X 


CTo/fb 


4.1232(2) 


4.1229(2) 


1.82773(5) 


4.1480(1) 


1.6063(2) 


(jLO/fb 


3.6930(5) 


3.6926(5) 


1.6484(1) 


3.7476(4) 


1.4355(4) 




-7.8 


-7.8 


-7.2 


-4.1 


-0.9 




-7.3 


-7.3 


-6.3 


-4.1 


-3.5 


Sqcd/% 


+24.9 


+24.9 


+24.6 


+24.9 


+25.1 


{Kqcd - 1)/% 


+39.5 


+39.5 


+38.1 


+38.2 


+40.0 




+0.3 


+0.3 


+0.0 


-0.0 


+1.0 




4.7884(5) 


4.7872(5) 


2.1332(1) 


4.9696(3) 


1.9566(4) 


_NLO / fu 
'^HAWK/ 


4.8635(5) 


4.8622(5) 


2.1616(1) 


5.0115(3) 


1.9706(4) 



Table 1: Born cross sections (Jq and a evaluated with NLO and LO PDFs, respectively, 
various relative corrections, and the resulting predictions for inclusive production in the various 
Higgs-strahlung processes (Mh = 120 GeV, Tevatron) according to ()3.8p and (j3.9p . All relative 
corrections 6 are given relative to ctq , and -K^qcd = ^qgd I ^^'^ ■ 



channel 


H/+zy; + X 


m-vi^x 


m+/- + X 


^vm + X 


^viln + X 


CTo/fb 


1.50846(7) 


0.66292(3) 


0.35349(2) 


0.74759(3) 


0.058236(9) 


CjLO/fb 


1.4183(2) 


0.60926(9) 


0.32845(5) 


0.69519(9) 


0.05417(3) 




-14.2 


-14.0 


-10.9 


-6.9 


-12.5 




-13.3 


-13.0 


-9.0 


-6.9 


-14.5 


<5qcd/% 


+9.5 


+9.4 


+9.8 


+9.8 


+6.8 


(i^QCD - 1)/% 


+16.5 


+19.1 


+18.1 


+18.1 


+14.9 




+1.3 


+1.5 


+0.0 


+0.0 


+12.5 




1.4522(4) 


0.6406(2) 


0.35329(7) 


0.7646(2) 


0.06043(6) 


"^HAWK/ 


1.4713(4) 


0.6488(2) 


0.35639(7) 


0.7697(2) 


0.06100(6) 



Table 2: Born cross sections, various corrections, and the resulting predictions for boosted pro- 
duction in the various Higgs-strahlung processes (Mh = 120 GeV, LHC at 7TeV). See Tabled] 
for details. 
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Figure 5: EW corrections with (left) and without (right) lepton-photon recombination for H/^z^; 
production at the LHC at 7 TeV and at 14 TeV and at the Tevatron as a function of the Higgs- 
boson mass. 



Figure [7] shows the differential cross sections with respect to the transverse momentum of 
the Higgs boson up to 300 GeV, along with the EW corrections. For the inclusive Higgs-boson 
sample, the EW corrections range between —5% and —15% and are largest for the W-mediated 
channels at large pt,h- The difference due to the treatment of final-state photons is again small, 
since the bulk of the correction is of weak origin and not due to final-state radiation. The same 
observation holds for large values of the transverse momentum in the boosted-Higgs analysis. 
However, close to the cut value of 200 GeV for pt,u, the corrections receive a sizeable negative 
contribution from final-state radiation effects, in particular for the H/+/^ channel, since we also 
require pt,f > 190 GeV. Final-state radiation can destroy the perfect balance between the Higgs- 
and the vector-boson transverse momentum at Born level. This balance is also influenced by 
initial-state radiation and would lead to a breakdown of the fixed-order calculation near the cut if 
the cuts on pt,h and Pt,w/z were equal. As mentioned above, the photon-induced contributions 
are large for the Hvi/Vi channel in the boosted analysis, but the overall contribution is small and 
phenomenologically unimportant. 

Figure El shows the differential cross section with respect to the transverse momentum of the 
vector boson, i.e. the sum of the transverse momenta of the leptons in the final state (or the 
transverse momentum of the neutrino in the Hz/; /Vi discussed at the end of Section [2]). 

The EW corrections are very similar to the ones for the pt,r distribution. Only the large 
contribution due to final-state radiation near the cut in the boosted analysis is absent because 
the Higgs boson is uncharged (possible final-state radiation of the Higgs decay products is not 
available in our analysis). However, near Pt,w/z = 200 GeV the perturbative uncertainty due 
to QCD initial-state radiation rises, since there is no tree-level distribution below the value of 
the pt,h cut. In Figure [U one observes a sharp drop of the cross section in the first bin of the 
distribution which is due to correspondingly large NLO QCD corrections. 

The differential distributions with respect to the rapidity of the Higgs boson and the trans- 
verse momentum of a positively charged lepton in the final state are displayed in Figures [9] and 
\T0\ respectively. The EW corrections for the rapidity distribution are almost completely flat and, 
hence, resemble the correction for the cross section using a given set of cuts. The transverse- 
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Figure 6: Relative correction from photon-induced processes for basic cuts (left) and the boosted- 
Higgs analysis (right) at the 7TeV LHC as a function of the Higgs-boson mass. 

momentum distribution shows corrections growing in absolute size with the energy of the typical 
events as observed before, close to or even beyond —20% at Pt,/+ = 300 GeV with recombination 
or for bare muons, respectively. 

Finally, Figure [TT] shows that for the H/^z^; channel the EW corrections hardly differ for the 
same observable measured at different hadron colliders. Only the photon-induced corrections 
are more sensitive to the specific collider settings. Exemplarily, the distribution with respect to 
the transverse momentum of the Higgs boson is shown but other distributions show the same 
qualitative behaviour. 

4 Conclusions 

In this paper we have discussed the impact of electroweak radiative corrections to the Higgs- 
strahlung processes off W/Z bosons, pp/pp — )• Yilvi / 1~ / vivi + X, at the Tevatron and the 
LHC. Compared to a previous calculation, we provide fully differential predictions. We include 
the leptonic decays of the W and Z bosons and support the full kinematical information on 
the Higgs boson and the decay leptons of the weak gauge bosons. These processes and the 
corresponding next-to-leading-order QCD and electroweak corrections are included in the new 
version 2.0 of the Monte Carlo program Hawk, which was initially designed for Higgs + 2 jet 
production via vector-boson fusion and WH/ZH production with hadronically decaying W/Z 
bosons. 

The electroweak corrections are of the order of —(5—10)% for total cross sections and, thus, 
larger than the uncertainty originating from parton distribution functions and QCD factorization 
and renormalization scales. Our detailed discussion of the electroweak corrections to transverse- 
momentum and rapidity distributions of the Higgs boson and outgoing charged leptons demon- 
strates that the corrections further increase in size in the regions of increasing transverse mo- 
menta. For instance, imposing the cut pt,h > 200 GeV at the LHC, which is part of the search 
strategy for highly-boosted Higgs bosons decaying into "fat jets" with bb substructure, drives the 
corrections to WH production to about —14% for Mh = 120 GeV. The electroweak corrections 
to the different distributions can already exceed —15% in the range = 200— 300 GeV and 
grow further when more energetic events are analyzed. 
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Figure 7: The pt,r distributions daf^[^ /dp^^n (top), NLO EW corrections for recombined and 
bare leptons (middle), and corrections due to photon-induced processes (bottom) for Higgs- 
strahlung with basic cuts at the Tevatron (left) as well as for boosted Higgs bosons at the 7 TeV 
LHC (right) for Mr = 120 GeV. 
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Figure 8: The pt,v distributions daf^[^ /dp^y (top), NLO EW corrections for recombined and 
bare leptons (middle), and corrections due to photon-induced processes (bottom) for Higgs- 
strahlung with basic cuts at the Tevatron (left) as well as for boosted Higgs bosons at the 7 TeV 
LHC (right) for Mr = 120 GeV. 
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Figure 9: The yu distributions daf^^ /dyu (top), NLO EW corrections for recombined and bare 
leptons (middle), and corrections due to photon-induced processes (bottom) for Higgs-strahlung 
with basic cuts at the Tevatron (left) as well as for boosted Higgs bosons at the 7TeV LHC 
(right) for Mr = 120 GeV. 
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Figure 10: The PT.i+ distributions daf^^ /dpT^ i+ (top), NLO EW corrections for recombined 
and bare leptons (middle), and corrections due to photon- induced processes (bottom) for Higgs- 
strahlung with basic cuts at the Tevatron (left) as well as for boosted Higgs bosons at the 7 TeV 
LHC (right) for Mr = 120 GeV. 
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Figure 11: The pT,H distributions daf^^/dpT^n (top), NLO EW corrections for recombined and 
bare leptons (middle) and corrections due to photon-induced processes (bottom) for for Hl'^i'i 
production and different cohider setups with basic cuts (left) and additional cuts for boosted- 
Higgs analyses (right) for Mh = 120 GeV. 
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Electroweak corrections, thus, represent an important ingredient in the theory predictions for 
Higgs-strahlung off W/Z bosons at the Tevatron and the LHC. 
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